[gm, gmpos] = max(g)
[g2m, g2mpos] = max(convol(g,g))
g2 = convol(g,g);
g2(40:60)
